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ABSTRACT 

We present OVRO Millimeter Array A = 3.4 — 1.2 mm dust continuum 
and spectral line observations of the accretion disk encircling the T Tauri star 
LkCa 15. The 1.2 mm dust continuum emission is resolved, and gives a minimum 
diameter of 190 AU and an inclination angle of 57 ± 5°. There is a noticeable, 
but at present poorly constrained, decrease in the continuum spectral slope with 
frequency that may result from the coupled processes of grain growth and dust 
settling. Imaging of the fairly intense emission from the lowest rotational tran- 
sitions of CO, ^^CO and HCO+ reveals a rotating disk substantially larger than 
that observed in the dust continuum. Emission extends to ~750 AU and the 
characteristic radius of the disk is determined to be ~425 AU (HWHM), based 
on model fits to the CO velocity field. The measured line ratios demonstrate 
that the emission from these species is optically thick, while that from C^'^O 
and H^^CO"*" is optically thin, or nearly so. The disk mass derived from the CO 
isotopologues with "typical" dense cloud abundances is still nearly two orders 
of magnitude less than that inferred from the dust emission, the most probable 
explanation being extensive molecular depletion in the cold, dense disk midplane. 
Thus, while CO, HCO+ and their isotopologues are excellent tracers of the disk 
velocity field, they are not reliable tracers of the disk mass. 
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N2H+ 1— >0 emission has also been detected which, along with HCO+, sets 
a lower limit to the fractional ionization of 10~*^ in the near-surface regions of 
protoplanetary disks. This first detection of N2H"'" in circumstellar disks has 
also made possible a determination of the N2/CO ratio (~ 2) that is at least an 
order of magnitude larger than those in the envelopes of young stellar objects 
and dense clouds. The large N2/CO ratio indicates that our observations probe 
disk layers in which CO is depleted but some N2 remains in the gas phase. Such 
differential depletion can lead to large variations in the fractional ionization with 
height in the outer reaches of circumstellar disks, and may help to explain the 
relative nitrogen deficiency observed in comets. 

Subject headings: stars: individual (LkCa 15) — stars: circumstellar matter — 
planetary systems: protoplanetary disks — radio lines: stars — ISM: molecules 

1. Introduction 

LkCa 15, located in the Taurus-Auriga complex at a distance of 140 pc, is rapidly 
becoming one of the best studied T Tauri stars at millimeter wavelengths. A K5 star with 
an optical/near- infrared luminosity of ~0.8 Lq (Elias 1978, Neuhaeuser et al. 1995, Simon 
et al. 2001), LkCa 15 is well isolated from dark clouds and lacks any detectable reflection 
nebulosity in HST WFPC2 images (Krist et al. 1997). Its substantial near- through far- 
infrared excess is consistent with the presence of a circumstellar disk (Strom et al. 1989) 
whose mass and radial extent have been determined by recent spectral energy distribution 
fits and millimeter-wave interferometric observations to be some 0.024 Mq and 650 AU, 
respectively (Chiang et al. 2001; Simon et al. 2001). 

Age estimates for LkCa 15 range from 3 Myr (Simon et al. 2001) to 12 Myr (Thi 
et al. 2001). These time scales are similar to those suggested for the disappearance of small 
dust grains in the inner regions of circumstellar disks (Skrutskie et al. 1990; Lada 1999; 
Robberto et al. 1999) and for the differentiation of meteoritic material in the early solar 
nebula (Alexander et al. 2001; Amelin et al. 2002). 

The large size and mass of the disk, combined with the age of LkCa 15, make this 
an important system for further study since it may represent an important transitional 
phase in which viscous disk spreading and dispersal competes with planetary formation 
processes. Further, the outer disk can be studied in detail with existing (sub) millimeter- 
wave telescopes. Single dish observations have already revealed readily detectable emission 
from several molecules (van Zadelhoff et al. 2001; Thi et al. 2001), while aperture synthesis 
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imaging has the potential to shed hght on the chemical and physical gradients in the disk. 

Our initial millimeter- wave observations of LkCa 15 have concentrated on species that 
are typically abundant in dense molecular clouds. Here we present Owens Valley Radio 
Observatory (OVRO) millimeter array observations of the disk encircling LkCa 15 in dust 
continuum emission and the lowest rotational transitions of various isotopologues of CO and 
HCO"'". An upper limit to the emission from the H^^CO"'" ion is noted. N2H"'" was also 
observed to constrain the fractional ionization and molecular depletion beyond 100 AU. Due 
to the limited angular extent of the disk, small scale gradients were largely smoothed out 
by the spatial resolution achieved, and the resultant line intensities should be considered 
as an "average" over the outer disk. Nevertheless, these interferometric observations suffer 
from substantially less beam-dilution than those with single telescopes; and this, combined 
with the inherent stability of cross-correlation measurements, provides a powerful means of 
studying the weak emission from species such as N2H+ and C^^O. While arrays under devel- 
opment will provide substantially improved sensitivity and spatial resolution, observations 
at the ^few arcsecond level form the first and necessary step to understanding depletion and 
ionization within the disks of T Tauri stars. Future publications will present more detailed 
analyses of these and other data acquired in an extensive millimeter-wave interferometric 
spectral line survey of LkCa 15 (Qi, Kessler & Blake 2003; Kessler, Qi & Blake 2003). In 
what follows, §2 includes a brief overview of the observations, while in §3 we discuss the 
results and the important effects of gas-phase depletion and fractional ionization on the 
physical and chemical structure of the disk around LkCa 15. A brief summary of our results 
and the role of future observations may be found in §4. 



2 . O bser vat ions 

All measurements were made between 1997 October and 2000 April using the OVRO 
MiUimeter Array at Big Pine, Cahfornia. The array consists of six 10.4 meter telescopes with 
an rms surface precision of 35 microns. The pointing accuracy is about 4", except for brief 
excursions at sunrise and sunset. Combinations of five array configurations were used to 
map LkCa 15, with antenna spacings ranging from 20 to 400 m east-west and from 25 to 440 
m north-south. Cryogenic SIS receivers on each telescope produced average single sideband 
system temperatures of 450, 1300, and 1200 K at the frequencies of the ^'^CO 1—^0, CO 1—^0, 
and CO 2^1 lines, respectively. The receivers were tuned for double sideband operation so 
that both the upper and lower sidebands could be used for molecular line studies. 

The receiver IF output is sent over fiber optic transmission lines to a pair of 1 GHz 
bandwidth analog correlators and a 512 lag digital correlator with four separately config- 
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urable modules. Spectral line and continuum measurements are made simultaneously and 

can be carried out in both the A = 3 and 1.3 mm atmospheric windows with a single local 
oscillator set up. During the ^^CO 1^0 and HCO+ 1— >0 observations of LkCa 15, the digital 
correlator was configured into two bands of 64 x 0.125 MHz channels to provide a velocity 
resolution of 0.17 and 0.21 km s~^, while the CO 2—^1 measurements were made using 256 
X 0.125 MHz channels for a velocity resolution of 0.16 km s~^. 

Calibration of the visibility phases and amplitudes was achieved with observations of 
the quasars 0528+134 and 0507+179, typically at intervals of 20-30 minutes. Measurements 
of Uranus, Neptune, 3C 273 and 3C 454.3 provided an absolute scale for the calibration 
of flux densities. All data were phase- and amplitude-calibrated using the MMA software 
package, developed specifically for OVRO (Scoville et al. 1993). Continuum and spectral 
line maps were generated and CLEANed using the NRAO AIPS package. The continuum 
level is not subtracted in obtaining the line maps because most of the transitions observed 
here are optically thick whereas the continuum is optically thin. Uncertainties in the fiuxes 
and source positions are estimated to be 20% and 0."5, respectively. The synthesized beams 
vary from 1" at 1.2 mm to 4" at 3.4 mm. 



2.1. Continuum Emission 

Figure 1 shows the spectral energy distribution of the continuum emission from LkCa 15. 
A noticeable, but at present poorly constrained, decrease in the continuTim spectral slope 
with frequency has been found, as discussed further in §3.1. The dust emission is spatially 
resolved, but only in the highest resolution OVRO 1.2 mm images. The (■u,f)-plane fits 
of a two-dimensional gaussian to the highest quality visibilities yield a deconvolved size of 
(l."35 ± 0."12) X (0."74 ± 0."06) at a position angle (PA) of 66 ± 5°, which corresponds to a 
minimum size of the dust disk of 189 ± 17 AU. Assuming the disk is circular, the measured 
ellipticity indicates an inchnation angle of 57 ±5° (0° being the inclination of a face-on disk). 
This inclination is somewhat larger than the value obtained by Duvert et al. (2000), 34°, 
but smaller than the value of 66° obtained by Kitamura et al. (2002) with their NMA 2 mm 
continuum observations, probably due to differences in the synthesized beams. Nevertheless, 
it is in good agreement with both the OVRO and Plateau de Bure measurements of the disk 
properties as traced by CO (Simon et al. 2001). 
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2.2. CO/HCO+ Observations 

Figure 2 shows the integrated intensity and intensity- weighted mean velocity maps of the 
CO 2— >1 emission toward LkCa 15. As shown in the left panel, the molecular gas emission 
appears much more extended than the dust. Solid grains no doubt exist at similar radii, but 
are undetectable due to the more rapid fall-off in dust emissivity (Sargent & Beckwith 1987; 
Lay et al. 1997). Conversely, the molecular emission from the innermost region of the disk 
(R<30-50 AU) is sufficiently beam diluted even at 1-2" resolution that no direct information 
on this size scale can be obtained without further improvements in sensitivity and imaging 
performance. The circumstellar gas emission is an effective tracer of the line-of-sight velocity 
component of the outer disk via Doppler shifts, as is shown by the intensity-weighted mean 
gas velocity map in Figure 2 (right panel). Molecular gas southwest (SW) of the star is blue- 
shifted, whereas that to the northeast (NE) is red-shifted, consistent with gas orbiting in a 
disk inchned with respect to the plane of the sky. With a beam size of ~2" (or 280 AU), the 
OVRO CO 2— >1 image of the LkCa 15 disk is well resolved. A gaussian fit to the velocity- 
integrated CO intensity map (Figure 2) results in a FWHM of (4."4 ± 0."4) x (3."6 ± 0."3) 
at a position angle of 68°, which corresponds to a radius of ^300 AU. 

At such moderately large inclinations shear broadening exceeds thermal broadening in 
the spectrum, and the local emissivity of optically thin lines becomes anisotropic (Horne & 
Marsh 1986). In this case, the separation of the two peaks in the CO 2— >^1 spectrum (2.6 km 
s"^, see Figure 3) becomes 2VfeepSini, where V^ep is the Doppler shift due to the Keplerian 
orbital motion, or indicating Ra^SlS AU. In Figure 4, observed maps on velocity 

channels of width 0.65 km s~^ are compared with synthetic maps generated by a detailed 
model of a disk in Keplerian rotation. Following Koerner (1995) (see also Mannings et al. 
1997), the model used a two-component gaussian emissivity distribution with half- maximum 
radii taken from the fit to the integrated intensity map in Figure 2. The inclination angle i, 
outer cut-off radius Rd, and relative amplitudes of the gaussian components were varied to 
match the structure in Figure 4. Best fits were found for a central mass 0.9 ± 0.2 Mq and 
disk inclination of idisk — 58° ± 10°. The outer gaussian component has a HWHM of 425 
AU and a cut-off radius of 750 AU. These values are in good agreement with optical and 
infrared photometry (Bouvier et al. 1995), PdBI CO observations (Simon et al. 2001) and 
with our resolved 1.2 mm continuum observations. As can be seen in Figure 4, synthetic 
maps produced from this model agree very well with the CO J=2-^l spectral line maps 
of LkCa 15. The size of the gas disk is thus well determined both from the image and the 
spectral data cube, and is much larger than the dust disk size as determined from continuum 
emission. In what follows, important disk physical properties are fixed at the values derived 
from the dust (disk mass) and CO (disk size and orientation) emission. 
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For robust measurements of the abundance and excitation of CO, it is necessary to 
derive observational estimates of line opacities. These can be obtained by examining a 
range of transitions and isotopologues. Visibilities from several configurations were therefore 
combined to examine the ^^CO, C^^O, HCO+, and N2H+ J=1^0 transitions; results are 
presented in Figure 5 and Figure 6. The velocity patterns observed in both ^^CO and HCO"'" 
are entirely consistent with those seen in CO, with red-shifted emission to the NE and blue- 
shifted emission to the SW. The C^^O and N2H+ 1— >0 data cubes have lower signal-to-noise 
values, and so no velocity analysis was carried out for these transitions. The integrated line 
intensities and the corresponding column densities for the above transitions were derived 
with the assumptions of local thermodynamic equilibrium (LTE) and optically thin emission 
at an excitation temperature of 30 K as presented in Table 1. Only upper limits have been 
obtained for Hi3cO+. 



3. Discussion 

3.1. Disk Structure and Grain Properties 

Two of the most important global properties of a disk are its mass and temperature 
distribution. At sufficiently short wavelengths, the entire disk is optically thick, and the 
density distribution does not affect the emission. At wavelengths longer than about 300 
yum, however, the thermal emission from circumstellar disks becomes optically thin. Thus, 
continuum observations at millimeter wavelengths provide an excellent way to measure disk 
masses directly via the equation: 

f r f a [ '^^T ^ n 2k < T > ^ ^ 

= / i,dn^ / s^T^n^ / -j^K, ^^^^ xw^ ' 

where M = J Edcx, F^, is flux density, Hi, is the mass opacity, S is the surface density and a 
is the surface area of the disk (Beckwith et al. 1990). The dust temperatures T are almost 
always sufficiently high (typically 50 K) to make the Ray leigh- Jeans assumption vahd, so the 
total emission is proportional to a product of the total mass and the average temperature 
< T >, weighted appropriately by its radial distribution. The major uncertainty in mass 
estimates arises from k^,. Theoretically, is expected to vary as a power of the frequency u: 

i^v = f^o ■ Following Beckwith et al. (1990), we adopt a fiducial value of 0.02 { 230a hz ) 
cm^g^^ and calculate the total gas mass assuming Mgas/Mdust — 100. The disk mass of 
LkCa 15 is calculated to be ~0.01 M0. 

More detailed constraints for the density and temperature structure can be derived 
by fitting the spectral energy distributions (SED) of young star/disk systems. Figure 1 
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shows the flux density as a function of frequency of the miUimeter wavelength emission from 
LkCa 15. Because the continuum and molecular lines were observed simultaneously, the 
frequencies of the continuum observations at 1.3 and 2.7 mm vary slightly depending on 
the local oscillator frequency required for the molecular line observations. In Figure 1, the 
spectral index, a, as defined by ^ z/", is 2.65, consistent with IRAM PdBl observations 
(Duvert et al. 2000). A separate analysis of the data in each atmospheric window, however, 
suggests that the slope of the 3 mm points alone {a — 3.51) is considerably steeper than 
that for the 1.3 mm fluxes {a — 2.34), although the error bars are large. 

In order to examine this issue further, wc compare our OVRO 3.4 - 1.2 mm continuum 
fluxes with those obtained at 870 /im and 7 mm by the JCMT and VLA, respectively. The 
JCMT SCUBA flux hes along a spectral index a of 2.65 with respect to the OVRO data, and 
so is consistent with the local slope derived from the OVRO 1.3 mm continuum emission 
measurements. No simple extrapolation of the sub-millimeter data can explain the VLA 
non-detection at 7 mm (D. Wilner, private communication), which requires a values similar 
to or even larger than those measured at 3 mm. The coincidence of the 1.3 and 3 mm data 
observed by Duvert et al. (2000) with the fits to our OVRO data indicate that the observed 
change in the spectral index at millimeter wavelengths is not likely to be due to errors in 
flux calibration. 

Several physical mechanisms, notably optical depth effects and grain growth, can con- 
tribute to the observed behavior. As mentioned above, /tj, is expected to vary as u^, where j3 
is the opacity index. Assuming optically thin emission in the Raylcigh- Jeans limit, the opac- 
ity index (3 {=a — 2) varies from 1.51 to 0.34 for the LkCa 15 disk, while f3 is approximately 
2 for dust in the diffuse interstellar medium. 

How reasonable is this assumption of optically thin continuum emission? Following 
Beckwith et al. (1990): 

where ,^ , ^ — r is the ratio of optically thick to optically thin emission and depends on 

V V (2— p)r 

the disk structure, = ^^^^^^2 is the average optical depth at frequency ly, and Mn is the 
mass of the disk. Given the dust disk radius of 190 AU and an Md of 0.01 M0 (derived 
above from F,^), the average optical depth at 230 GHz is ^0.03. Thus, for any reasonable 
disk density and temperature distribution, the optical depth correction does not appear to 
be sufficient to explain the observed change in the LkCa 15 spectral index with a single /3. 

One means to solve this dilemma is for the mass opacity coefficient K,,y to decrease by 
nearly an order of magnitude between 1.2 and 3.4 mm. Recent detailed models of irradiated 
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T Tauri disks, including dust grain growth (D'Alessio et al. 2001), conclude that when grains 

grow from yum to mm sizes, the mass opacity coefficient changes significantly between 1.3 
and 3 mm. Indeed, observations of circumstellar disks (Beckwith ct al. 1990; Mannings & 
Emerson 1994; Wilner et al. 1996) indicate typical values of beta around 0.3, similar to 
that expected from mm-sized pebbles (0 < /9 < 2), as opposed to the steeper frequency 
dependence of Ky {j3 = 2) found in the general interstellar medium (Hildebrand 1983). 

As indicated by Chiang et al. 2001, the disk mass and grain size distribution cannot be 
uniquely constrained by the continuum SED alone when only unresolved fluxes are available. 
Additionally, D'Alessio et al. (2001) flnd that for a given total dust mass, grain growth is 
accompanied by a decrease in the vertical height of the disk surface. Fits of the SED for 
LkCa 15 (Chiang et al. 2001) indicate that the ratio of the height of the disk photosphere 
(H) to the vertical gas scale height (h) is small (1.0), which can be interpreted as evidence 
that the dust has begun to settle vertically toward the midplane. The H/h parameter is a 
rather robust one within the confines of the model, for as H/h decreases from 4 to 1 there is 
an corresponding decrease in the overall level of infrared excess at A <100 //m since the disk 
intercepts less stellar radiation. Dust settling can help explain the relative paucity of large 
disks to which the powerful technique of scattered light coronography can be applied and 
specifically the lack of any detectable refiection nebulosity in HST WFPC2GT0 observations 
of LkCa 15 (Krist et al. 1997). Thus, the most hkely explanation for the observed change in 
the spectral index from 1 mm to 3 mm is a combination of grain growth and dust settling. 



3.2. Radiative Transfer and Moleculcir Line Emission 

The observed line intensities in circumstellar disks are a complex function of the physical 
structure of the disk and of the line/ continuum optical depth. If we assume that circumstellar 
disk abundances are similar to those found in dark clouds, the emission lines from parent 
isotopologues of common species such as CO, HCO+, HCN, etc., should be highly optically 
thick. Our observed ratios of ^^CO/^^CO emission reveal that ^^CO is indeed optically 
thick, but say httle about ^^CO. Several LkCa 15 transits were therefore devoted to a local 
oscillator setting containing the C^^O J=l transition, the results of which are presented 
in Figure 5. The observed beam-matched flux density ratio is 5, only somewhat lower than 
the isotope ratio [^^C]/[^^0] of 8.3. The ^^CO J=l transition is therefore only somewhat 
optically thick, and that of C^^O is optically thin, and should therefore provide a more reliable 
estimate of the disk mass, provided it remains in the gas phase, van Zadclhoff ct al. (2001) 
indicate that, for standard HCO^ abundances, the 1— >0 line is again optically thick in the 
outer layers, whereas that of H^^CO+ is close to being optically thin throughout the disk. 
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Further, the low-J HCO+ lines are predicted to probe densities of order 10^ — 10^ cm~^, 
below the critical density of the 4— transition (so that LTE cannot be safely assumed). 
For H^^CO"'", the populations will be closer to thermal equilibrium because its emission arises 
primarily from regions with densities of 10^ — 10^ cm~^ (van Zadelhoff et al. 2001). Parent 
isotopologue lines are therefore likely to be optically thick for many species in disks, and 
so whenever possible the emission from isotopically substituted species should be utilized in 
deriving relative molecular abundances. 

The observations described above indicate the need for detailed calculation of disk ra- 
diative transfer. Furthermore, the complex geometry within the disk demands that flexible 
and efficient radiative transfer techniques be utilized, and of the many techniques available 
the Monte Carlo approach has shown great promise (Choi et al. 1995). Since Monte Carlo 
techniques utilize integration paths chosen at random, any number of coordinate systems 
can be used. The high optical depth of many transitions in circumstellar disks generates 
very short photon mean free paths, however, and so unless aggressive acceleration schemes 
are employed, the Monte Carlo approach can be very slow. Computational instabilities and 
poor convergence can result, but efficient one-dimensional (ID) and two-dimensional (2D) 
Monte Carlo treatments of the line radiative transfer based on the two-layer passive disk 
models outlined above have now begun to appear (Hogerheijde & van der Tak 2000; van 
Zadelhoff et al. 2001). In all such approaches, the intensity of each fine is calculated by 
solving the radiative transfer equation. 

2D implenetations are necessary to quantitatively treat inclined disks. For this reason 
we have used an accelerated 2D Monte Carlo model (Hogerheijde & van der Tak 2000) 
to examine the radiative transfer and molecular excitation in the LkCa 15 disk, taking 
both coUisional and radiative processes into account. This model produces a simulated 
observation of each transition as imaged by a telescope with infinite resolution for a disk of 
a given size, inclination and temperature distribution. The MIRIAD function UVMODEL 
and the observed visibility data set were then used to sample this model at the observed 
(u, v) spacings and the model data set was processed in a manner identical to that of the 
OVRO data, thus allowing a direct comparison of the two. For these model calculations, the 
temperature and hydrogen density as a function of radius and height were acquired from a 
model similar to those presented in D'Alessio et al. (2001), but calculated specifically for 
the stellar parameters of LkCa 15 such that the overall SED of the star-|-disk is reproduced. 
In this model, we assume spherical, compact dust grains with a power law size distribution 
with exponent 3.5, a maximum grain size of 1 mm and grain abundances from Pollack et al. 
(1994). The accretion rate is assumed to be 1.0 x 10~* Mq yr~^, with the inner radius of the 
disk set at 5 AU and the outer radius at 425 AU, which is consistent with that derived above 
for LkCa 15. The inclination (58 degrees) and turbulent velocity width (0.1 km s~^) are found 
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by fitting the observed CO 2— >1 line profile with the model spectrum, using Nco—^^'^^h ^ 
as shown in Figure 3. For each transition observed, a suite of models with the parameters 
described above were created, but with the molecular abundances relative to hydrogen varied 
until a match with the observed integrated line intensity was obtained. Integrated intensity 
is a useful quantity for comparison because the column density is approximately constant as 
a function of radius (Willacy & Langer 2000; Aikawa & Herbst 1999) at the large radii (> 
50 AU) to which we are sensitive with the OVRO millimeter array. 

The resulting column densities are shown in Table 1. The column densities calculated 
using this non-LTE radiative transfer model arc 1-2 orders of magnitude larger than those 
calculated in the Willacy & Langer (2000) models of disk chemistry, but consistent with 
those calculated by Aikawa et al. (2002) - except that N2H+ is an order of magnitude larger 
than the model prediction. As indicated in the last column of Table 1, for CO 2— >^1, ^^CO 
1— >-0, HCO+ 1— >0 and N2H+ 1— >0 emission the column densities are smaller when calculated 
using the LTE assumption. This is due to the fact that the emission from these transitions is 
optically thick (see discussion above). The LTE calculations assume that the gas is optically 
thin and that the entire disk is being probed by observation of this line emission and thus 
result in an underestimation of the total amount of CO, HCO"'" and N2H+ present. In the 
case of the C^^O 1^0 and H^"^CO+ 1— >0 emission, which are believed to be optically thin, 
and for the conditions in the near-surface regions of disks beyond 50 AU, the column densities 
are overestimated when thermal equilibrium is assumed. This supports the conclusion that 
conditions in the emitting region are not at LTE, and that both radiative and coUisional 
processes play a part in the excitation. The population of higher energy rotational states 
is thus much lower than would be the case if collisions were dominant; and the assumption 
of thermal equilibrium therefore results in an overestimation of the population of the upper 
states and thus in the total column density, which at LTE is directly proportional to the 
density of the upper state. 

Both the results presented above and observations of high- J transitions of CO, HCO"'", 
and HCN toward LkCa 15 by van Zadelhoff el al. (2001) suggest that the molecular emission 
from the major isotopologues is optically thick and that the inferred densities in the region 
are not sufficient to thermalize all transitions. The line ratios presented in van Zadelhoff el 
al. (2001) further indicate temperatures and densities typically associated with warm layers 
lying above the disk midplane. Additionally, if we ignore depletion of gas phase species 
onto grains, and assume cosmic fractional abundances for CO (10~^-10~^ Nh) and hydrogen 



^This model assumes that CO and H2 have similar radial and vertical distributions. Chemical models 
Aikawa ct al. (2002)) indicate that this is not strictly the case, and variations in the vertical distribution 
of CO will be explored in future work. 
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column densities of 10^^ cm~^, the resulting CO column density is an order of magnitude 

greater than that observed here (10^^ cm~^). Thus, for any disk in which the gas and dust are 
well mixed and there is no depletion, even species such as C^^O are predicted to be optically 
thick at molecular cloud abundances. The low column densities derived here and the high- 
J/low-J line ratios therefore indicate that the disk midplane (heights below several tens AU) 
is not being probed by the observed emission and that the column densities presented in 
Table 1 are only applicable to the warm layers of the disk. 



3.3. Ionization Balance 

The theory of the ionization balance in molecular clouds and the relevant ion-molecule 
chemistry has been widely used for determinations of the electron abundances in dense clouds 
(Langer 1985). Theoretical chemical models indicate that ion- molecule reactions are also 
efficient in the outer regions of protoplanctary disks (Aikawa et al. 1998). Photoionization, 
cosmic ray ionization and, possibly. X-ray ionization of H2 rapidly lead to the formation of 
H3 , the pivotal species in ion- molecule models of molecular cloud chemistry. The low proton 
affinity of H2 ensures that H3 , which has no rotational spectrum, will undergo proton transfer 
reactions with abundant neutral species such as CO and N2. In disks, gaseous CO and N2 
exist in regions well removed from the midplane, where the temperature warms enough for 
sublimation of these molecules (T ^ 20-30 K). The isoclcctronic ions HCO+ and N2H"'" were 
observed here and can therefore be used to estimate the fractional ionization in the warm 
layer of the LkCa 15 disk. 

The volume density in the emitting region can be constrained from line ratios obtained 
from the results presented here and previous single dish observations of molecules with 
high dipole moments (e.g. HCO^, HCN). An LVG analysis for LkCa 15, performed by 
van Zadelhoff et al. (2001), indicates that the hues arise from regions with nH2~10^-10^ 
cm~^. We use 10^ cm~^ as the number density of H2 in the calculations below, which is also 
the number density at a height of ~150 AU and a radius of 350 AU in the D'Alessio disk 
model. The mean kinetic temperatures are less well constrained, but the analysis of multiple 
transitions of CO (van Zadelhoff ct al. 2001) suggests that Tfci„~20-40 K for LkCa 15. 
We therefore adopt a single kinetic temperature of 30 K. Also, for the reasons discussed 
in §3.2, we use the column densities resulting from the non-LTE radiative transfer model 
summarized in Table 1 for all observed species. Because H2 cannot be observed directly at 
millimeter wavelengths, the column density is highly model dependent. Typical values vary 
from 10^^ to 10^^ cm"^ for R~300 AU, the distance corresponding to the typical beam size 
of our observations, van Zadelhoff et al. (2001) investigate the level of molecular depletion 
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by statistical equilibrium calculations as compared to the gas column densities predicted by 

dust continuum observations, and find depletion by a factor of 3-40 for CO in the disk of 
LkCa 15. Wc therefore adopt 10~^ as the fractional abundance of CO, by assuming depletion 
by a factor of 10 from the nominal value of /(CO) in molecular clouds. In our calculations 
we also assume that the warm layers in the disk are in local (chemical) equilibrium and that 
the abundances of CO, HCO"*" and N2H+ peak within these layers (Aikawa et al. 2002), 
although the vertical distributions of these species are not identical [cf. Figure 4 of van 
Zadelhoff et al. (2003) and Figures 4 and 5 of Willacy & Langer (2000)]. 

To simplify the calculations presented here, wc first consider only the ionization balance 
determined by HCO+, H^, N2H+ and electrons in steady state, for which 

ki[H+][CO] = ke{HCO+)[HCO+][e] - k2[N2H+][CO] (3) 



keiH^mm + ksmm] + k,[Ht][co] = (4) 

h[N2][H+] = k2[N2H+][CO] + ke{N2H+)[N2H+][e] (5) 

where ke{H+) = 1.15 x 10"^ (^0)"°''' cmV^, ke{HCO+) = 2.0x 10"^ (^0)"°'^' ^m^ 
s"-*^ and ke{N2H'^) — 1.7 x 10^^ (^) cm^ s"-*^ are the dissociative recombination rate 
coefficients of H+, HCO+ and N2H+, ki = 6.56x 10-^° (^)~°'^ cm^ s-^ and ka = 8.8 x 

10"^'' cm^ s~^ are the reaction rate coefficients of CO with and N2H"'", respectively, ka = 
1.8 X 10~^ is the rate coefficient for the reaction N2 + H^ N2H^ + H2 (all rate coefficients 
have been adopted from the UMIST database for astrochemistry (Millar et al. 1997)), and 
the bracket [ ] represents the fractional abundance of the species and C, the ionization rate 
due to cosmic rays. X-rays, and UV radiation. 

A lower limit to the electron fraction follows directly from the observed HCO^ and 
N2H^ column densities, but can also be obtained more "locally" by noting that the right 
side of equation 3 must be larger than 0, or: 

[e] k2[N2H+] 

[CO] - ke{HCO+)[HCO+] ^ ^ ' 

The lower limit of the fractional ionization [e] is therefore ^10~^, given [CO]~10~^, in the 
near-surface region of the disk. This value is fairly high compared to that of cold, dark clouds 
at riHa^lO^ cm~^ where for example, Caselli et al. (2002) obtained fractional ionization of 
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10~^ in the center of L1544, and we note that the ionization fraction in LkCa 15 may be 
enhanced relative to other T Tauri stars. The ROSAT All-Sky Survey of young stars in 
Taurus (Neuhaeuscr ct al. 1995) shows that the probability of the existence of an X-ray 
source for LkCa 15 is 98% and the upper limit for its X-ray luminosity is 4.1 x 10^^ ergs^^. 
Ionization from the resulting high-energy radiation may have important effects on the disk 
(Feigelson & Montmerle 1999). The high HCO"^ abundance and CN/HCN ratio observed 
(Qi, Kessler & Blake 2003) indicate that LkCa 15 could possess similar chemistry to the 
intense X-ray source TW Hya, another Sun-like star with strong molecular emission from its 
attendant disk (Kastner et al. 1997). 

Understanding the ionization rate near the disk surface is important because the electron 
fraction scales linearly with the ionization rate. A full 3D calculation of X-ray transport and 
ionization in axially asymmetric disks has been developed by Igea & Glassgold (1999) using 
a Monte Carlo approach. They calculate the ionization rate at it! = 1 AU in a minimum 
mass solar nebula disk model with L^, = 10^^ erg s~^, which can be scaled to other X- 
ray luminosities and disk radii and is found to be relatively insensitive to changes in disk 
parameters. Scaling their values at i? = 300 AU with N(H2)~10^^ cm~^, we calculate an 
X-ray ionization rate of ~10^^^ s^^ near the surface of the LkCa 15 disk. This is nearly a 
factor of ten larger than the galactic cosmic ray ionization rate, and shows that X-rays can 
play an important role in the disk ionization balance. The enhanced ionization can drive 
additional chemical activity, and the inclusion of X-rays within the context of dust settling 
models developed for LkCa 15 may well account for much higher abundances of molecules 
observed toward this source as compared with the DM Tau disk, for which "standard" disk 
chemistry models produce good agreement with observations (Aikawa et al. 2002). 

By considering the balance of all the molecular ions present in the disk, even tighter 
constraints can be placed on the electron fraction. We consider two conditions with differing 
abundances of metal ions, as discussed in Langer (1985). When the abundances of metal 
ions arc low, Hjj', HCO^ and N2H"'" become the major cations within the disk. Because 
HCO+ and N2H+ have fractional abundances that are only iJ10~^° while the lower limit to 
the electron fraction is closer to 10~^, dominates. In this case, the electron abundance 
can be calculated simply from [e] ~ ( ^ n{H2) ^ ^ ^'^ ^ (Langer 1985), which yields 
values similar to that derived above. However, the high resulting abundance of drives 
rapid protonation of CO, that is h[H+][CO] > ke{HCO+)[HCO+][e], and makes HCO+ 
balance (equation 3) impossible within the confines of the simple model. 

It is therefore reasonable to examine the contribution of metal ions to the fractional 
ionization because they are destroyed very slowly and can comprise a significant fraction 
of the ion abundance even though their elemental abundances are small. Promang et al. 
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(2002) calculate the ionization fraction of protostellar disks, taking into account vertical 
temperature structure and the possible presence of trace metal atoms, and find that a tiny 
fraction of the cosmically available metals can dramatically affect the ionization balance. 
Considering the model above, but now with [H^] + [HCO^] + [A^2-f^^]^[e], we can solve 
for the 3 unknown variables [A^2], [H^] and [e] in order to balance equations 3,4,5. Using 
equations 4 and 5, removing the term A;3[iV2][i?^], and inserting the value of Hg derived from 
the equation 3, the equation for [e] becomes: 

ke{HCO+)ke{H+)[HCO+][e]' + B[e] - -^h[CO] = (7) 

where B = -k2keiH+)[CO][N2H+] + kih{HCO+)[CO][HCO+] + kike{N2H+)[CO][N2H+]. 
This equation can be solved using the observed column densities of CO, HCO^, and N2H"'" 
and an assumed density of H2. The N2/CO ratio can also be derived from equations 5 
and 7 under such conditions, and Figure 7 presents the solutions for a range of H2 column 
densities. Not surprisingly, the electron fraction is calculated to be larger than the metal 
free case (~10~^), which, as we stress above, can only be applied for the physical conditions 
prevailing near the warm surface layer(s) of the disk. 

Transport processes in the outer disk likely depend critically on the degree of ionization. 
Efficient coupling between the neutral material and magnetic fields is thought to be provided 
by fractional ionizations of ^ 10~^ (Feigelson & Montmerle 1999). As we have outhned above, 
photons. X-rays and cosmic rays can be important sources of ionization near the disk surface. 
The photon penetration depth is very short compared to the disk thickness, however, and 
even cosmic rays are exponentially attenuated for columns >100 g cm~^. The electron 
fraction can therefore be expected to drop toward the midplane due to the increasing density, 
suggesting the possibility of layered accretion (Gammie 1996), with a turbulent surface layer. 
However, due to the depletion of volatiles such as CO, the disk mid-plane fractional ionization 
may be larger than expected from a direct extrapolation of the electron abundance inferred 
for the surface layer. For very high depletions, electron recombination begins to control the 
fractional abundance of H3 , resulting in a higher than expected electron (and H3 ) abundance. 
The details of the vertical and radial ionization fraction are pivotal to magnetohydro dynamic 
(MHD) models of YSO accretion disk viscous transport mechanisms (Glassgold et al. 2000). 
Balbus & Hawley (1991), for example, propose that differential rotation drives instabilities 
in the perpendicular magnetic field of accretion disks. Such magnctorotational instabilities 
then create a turbulent disk if they are well coupled to the neutral disk material. As shown 
above, observations suggest that the combination of photons. X-rays, and cosmic rays can 
produce the ionization needed to couple the magnetic field to the neutral disk material near 
the disk surface. Chemically, transport or mixing processes between the surface layer and 
the disk interior would be necessary to bring freshly made molecules such as CN, HCN 
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or even deuterated species back to the interior of the disk. In this scenario, the extent 
of such processes, which are dependent on the ionization fraction, would greatly effect the 
composition of early solar system objects such as comets. 

3.4. CO Depletion and the N2/CO Ratio 

If CO remains primarily in the gas phase, it is sufficiently chemically stable that its 
abundance will be httle altered from dark cloud values. CO would therefore be a good tracer 
of the disk gas mass, but depiction onto grain mantles in the cold disk is likely to substantially 
reduce the gas phase column density of CO and many other species. Thus, estimates of the 
degree of depletion of volatile molecules such as CO and N2 onto grains will greatly improve 
our understanding of disk chemistry. For CO, depletion factors can be measured directly 
from infrared spectroscopy of edge on disks (Boogert et al. 2002), or more generally by using 
observations of the dust continuum along with several isotopically-substituted molecules 
(e.g. ^^CO, ^^CO, C^^O) to distinguish depletion from optical depth effects. Prom the column 
density of CO (1.7 x 10^^ cm~^) obtained with the detailed 2D Monte Carlo radiative transfer 
model (Table 1), a disk gas mass of only ~10~^ M0 is derived, a value in agreement with that 
derived from single dish observations of ^^CO higher-J transitions (Thi et al. 2001) using 
the JCMT and CSO but nearly two orders of magnitude less than the mass derived from 
dust emission - far more than can be accounted for by errors either in the dust mass opacity 
coefficient or in the spectral line radiative transfer. Depletion thus appears to be a dominant 
effect and while CO and its isotopomers are therefore excellent tracers of the disk velocity 
field, it is clear that they are not reliable tracers of the disk mass. In the models of Aikawa 
and co-workers (Aikawa et al. 1996, 1999), CO depletes onto grains at radii greater than 200- 
300 AU. Inside this radius, the temperature warms to the subhmation temperature of CO (T 
~ 20 K) and the gas phase column density rises substantially. If thermal sublimation is the 
main desorption mechanism, the depletion patterns in disks should be strongly correlated 
with volatility. We have observed strong emission from non-volatile species such as HCN 
and CN toward LkCa 15 with spatial distributions that are not well correlated with CO 
(Qi, Kessler & Blake 2003). Additional mechanisms are therefore likely to be important in 
maintaining the observed gas-phase abundances. 

As a closed shell homonuclear diatomic molecule, N2 has no electric dipole allowed 
rotational transitions and so cannot be observed directly at millimeter wavelengths. However, 
from our N2H"'" observations and the model presented above, we can estimate the N2/CO 
ratio, and find it to be ~2 (Figure 7). This is much higher than the ratios inferred for typical 
dark clouds at moderate densities, which are ~0. 01-0.1 (Womack et al. 1992; Bergin et al. 
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2001). In the dark cloud IC5136, Bergin et al. (2001) do find evidence for the presence of 
differential gas-phase depletion and altered CO/N2 ratios in the densest portions, where CO 
exhibits a significant abundance reduction and N2H+ is relatively undepleted - although in 
their models the abundance of N2 is still less than that of CO. 

Within disks, Aikawa & Herbst (1999) investigated the two-dimensional chemical struc- 
ture within the so-called Kyoto minimum mass solar nebula model (Hayashi 1981), and pre- 
dicted much lower abundances of N2 than CO at all vertical scales. Because they adopted 
an artificially low sticking probability to reproduce the spectra of CO available at that 
time, without specifying any non-thermal dcsorption processes or modifying the tempera- 
ture distribution in the Kyoto model, such an approach might not account for the known 
differential desorption of CO and N2. Further, the vertical temperature gradients of actual 
protoplanetary disks should be more complex than those in the Koyoto model (Chiang & 
Goldreich 1997; D'Alessio et al. 1997). Willacy & Langer (2000) investigated the molecules 
in the super-heated layer of the Chiang and Goldreich (1997) disk models, and found that 
molecules in this layer are destroyed by the harsh ultraviolet radiation from the star, al- 
though in the warm upper layers of the D'Alessio et al. models UV shielding is sufficient to 
prevent destruction of most molecules (Aikawa et al. 2002). We therefore believe that N2 
and its chemical product N2H+ trace dense gas in the disk where CO is largely depleted but 
N2 is not, and that this region is protected by an overlying warm layer that protects N2 from 
ultraviolet radiation, in order to produce the high N2/CO ratio observed. 

Deriving absolute fractional abundances is more difficult, but from the N2/CO ratio 
and our Monte Carlo radiative transfer calculations we estimate a fractional abundance for 
N2 of ~ 2 X 10~^. N2 would therefore represent a substantial fraction of the cosmically 
available nitrogen, this plus the gas phase N2/CO ratio inferred for LkCa 15 provides a 
direct observational connection to comets. The abundance of N2 in comets is an important 
guide to their volatile content because it is trapped and released by amorphous ice in a 
manner similar to that of Ar (Bar- Nun et al. 1988). This process has implications for 
understanding the role of comets for delivery of volatiles and noble gases to the terrestrial 
planets. Wyckoff et al. (1991) have reviewed the total nitrogen abundance in comet Halley 
and argued that it is depleted by 2-6 fold relative to that of the sun. The laboratory 
experiments performed by Bar-Nun et al. (1988) showed that CO is trapped 20 times more 
efficiently than N2 in amorphous ice formed at 50 K and Owen & Bar-Nun (1995) predicted 
that in icy planetesimals forming in the solar nebula at about 50 K, N2/CO ^ 0.06 in the 
gases trapped in the ice provided N2/CO 1 in the nebula itself, a value quite consistent 
with our data. Wyckoff & Theobald (1989) reported N2/CO ~ 2 x 10~^ in Hallcy's comet, 
whereas Lutz et al. (1993) found N2/CO ~ 2 x lO'^ in Comet Bradfield, and N2/CO - 
3 X 10~^ in Comet Kohoutek. The observations are actually of the two molecular ions (CO"*" 
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and NJ) and the conversion from the ion abundances to those of the neutrals is dependent on 
poorly constrained photodestruction branching ratios. Wyckoff & Theobald (1989) argued 
that the ion ratio must be multiplied by a factor of 2, while Lutz et al. (1993) find no factor 
necessary. Large amounts of gaseous N2 may thus be required in the outer solar nebula, as 
we find in the disk of LkCa 15, in order to explain the N2/CO ratios found in comets. This 
can be confirmed by measurements of the N2/CO ratio in dynamically new comets, which 
should have higher values than those in short-period comets (Owen & Bar- Nun 1995). 

4. Summciry 

We have detected molecular gas and continuum emission from the disk encircling the T 
Tauri star LkCa 15 at frequencies between 85 and 260 GHz. The 1.2 mm dust continuum 
emission is resolved in our l."35 x 0."74 beam with a minimum diameter of 190 AU and 
an inclination angle of ~57°. A noticeable decrease in the continuum spectral slope with 
frequency may result from the combination of grain growth and dust settling. Sensitive ob- 
servations at longer wavelengths and submillimeter continuum images are needed to confirm 
the trend. Based on Keplerian rotation model fits to the CO velocity field, the gas emission 
extends to ~750 AU while the characteristic radius of the disk is determined to be closer to 
425 AU from Gaussian fits to the GO 2^1 spectral data cube. 

We have also presented molecular images of various isotopologues of GO and the HGO+ 
and N2H"'" ions in the disk. Our observations show that most of the molecular emission from 
the dominant species of common molecules is optically thick. Different molecules and even 
different transitions of one molecule will therefore probe different layers of the disk; most of 
the emitting region detected is near the disk surface. Detailed radiative transfer modelling 
indicates that LTE is not appropriate in these regions. With this 2D model and the column 
densities from the detection of the HGO"*" and N2H+ ions, we derive a lower limit to the 
fractional ionization of 10~^ and an N2/GO ratio of 2, which reveals a region where GO is 
heavily depleted but N2 is not. This zone must be well protected by the efficient shielding of 
ultraviolet radiation by an overlying disk layer. Additional discussions concerning the UV 
field and its impact on disk chemistry will be presented in forthcoming papers. 
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Table 1. Observed molecular intensities and column densities toward LkCa 15 



Transition 


Beam 


/ Tdv 


N(30K,LTE) 


N(Model) 


Model Ratios 




(arcsec) 


(K km s"^) 


(cm~^) 


(cm~^) 


N(Model)/N(30K,LTE) 


CO 2^1 


1.81x1.49 


12.5 


7.28(15) 


1.68(18) 


230 


i^CO 1^0 


2.74x2.41 


6.39 


1.11(16) 


3.04(16) 


2.74 


C^^O 1-^0 


5.05x3.59 


1.90 


3.31(15) 


1.40(15) 


0.42 


HCO+ 1^0 


7.74x5.64 


3.30 


9.25(12) 


2.31(13) 


1.79 


H13C0+ 1^0 


6.30x4.41 


<0.88 


<2.60(12) 


<1. 12(12) 


0.43 


N2H+ 1^0 


3.51x3.17 


3.83 


1.71(13) 


3.07(13) 


1.80 
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Fig. 1. — Measured millimetcr/submillimeter continuum emission fluxes toward LkCa 15. 
The 850 point was acquired with SCUBA, whose image of LkCa 15 is consistent with 
a point source (G. Sandell, private communication). The two squares depict the 3 and 
1.3 mm continuum fluxes measured with the Plateau de Bure interferometer (Duvert et al. 
2000). The spectral slope as determined by all of the available data is shown by the solid 
line, individual spectral slope fits to the 3 mm and 1.3 mm data are depicted with the 
two dashed lines. 7 mm VLA observations (D. Wilner, private communication) detect no 
continuum emission to an rms of 0.29 mJy. Even the 3a upper limit is inconsistent with an 
extrapolation of the spectral index derived from the millimeter /submillimeter fluxes. 
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Fig. 2. — CO 2^1 emission toward LkCa 15. The integrated intensity is shown in contours 
at left, overlayed on the highest resolution continuum image in gray scale. Contours start at 
30% of the peak value of 2.66 Jy km s~^ and increase 10% thereafter. The intensity-weighted 
mean velocity field over the observed emission region is shown at right {vlsr units in m s^^). 
The synthesized beam is shown at lower left. 
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Fig. 3. — The disk averaged spectrum of CO 2— >1 (solid histogram), compared with the 
model CO emission for a 430 AU radius disk inchned at 58° (dotted line). Details of the 
calculation may be found in §3.2. 
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Fig. 4. — Spectral-line maps of the CO 2^1 emission toward LkCa 15, obtained with 
OVRO in steps of 0.65 km s~^, shown adjacent to simulations of the emission predicted by 
a kinematic model of a disk in Keplerian rotation with parameters as described in the text. 
The velocities are refered to the stellar (systemic) velocity of Visr = 6.3 km s~^. 



-26- 



LkCa T5-: ^^CO(I-O) 




22 15 06 



a (1950) 



22 15 04 



4^36'^19!'2 




18M 18^0 
a (1950) 



20 



o 

LO 



22 15 04 



1 1 1 1 

LkCa 15: C^^0(1- 

.'"1 \ 

- ^ 


3) 


1,1, 





16 



-J 2 - 



o 
in 
en 



•^8 



22 15 04 - 



4^36'^18^8 18^4 18^0 

a (1950) 

.UCa lk/NNH + (1 -0) 





4^36'^18^8 



18M IS^'O 17^6 

a (1950) 



Fig. 5. — Top: The integrated intensity maps of ^'^CO (left, contours start at 0.16 Jy beam"^ 
km and are spaced by 0.08 Jy beam~^ km s^"*^) and C^'^^O (right, contours start at 0.1 Jy 
beam~^ km s~^ and are spaced by 0.05 Jy beam^^ km s^^) 1^0 emission. Bottom: The 
integrated intensity maps of HCO+ (left, contours start at 0.16 Jy bearn"^ km s~^ and are 
spaced by 0.08 Jy beam""^ km s~^) and N2H+ (right, contours start at 0.1 Jy beam~^ km 
s~^ and are spaced by 0.05 Jy beam~^ km s~^) 1 — >^ emission. All steps correspond to the 
rms noise level in each map. 
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Fig. 6. — Disk averaged spectra of the molecules in Figure 5, obtained over boxes of 6"x7" 
for i^CO, 8"x5" for C^^O, 6"x7" for HCO+, and 6"x6" for N2H+. 
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Fig. 7. — Fractional abundances {f{X) = N{X)/N{H2) of N2 and electrons versus the H2 
column density for the simple fractional ionization model discussed in §3.3. 



